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Globular clusters in modified Newtonian dynamics: 
velocity-dispersion profiles from self-consistent models 
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ABSTRACT 

We test the modified Newtonian dynamics (MOND) theory with the velocity- 
dispersion profiles of Galactic globular clusters populating the outermost region of 
the Milky Way halo, where the Galactic acceleration is lower than the characteris- 
tic MOND acceleration a . For this purpose, we constructed self-consistent, spherical 
models of stellar systems in MOND, which are the analogues of the Newtonian King 
models. The models are spatially limited, reproduce well the surface-brightness pro- 
files of globular clusters, and have velocity-dispersion profiles that differ remarkably in 
shape from the corresponding Newtonian models. We present dynamical models of six 
globular clusters, which can be used to efficiently test MOND with the available ob- 
serving facilities. A comparison with recent spectroscopic data obtained for NGC2419 
suggests that the kinematics of this cluster might be hard to explain in MOND. 

Key words: gravitation - stellar dynamics - methods: analytical - stars: kinematics 
- globular clusters: general. 
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1 INTRODUCTION 

Modified Newtonian dynamics (MOND; Milgrom 1983) rep- 
resents one of the most popular alternative to the commonly 
accepted dark-matter paradigm. This theory was put for- 
ward in the early eighties as a way to explain the observed 
flat rotation curves of disk galaxies without the need of hid- 
den matter. In Bekenstein & Milgrom's (1984) formulation 
of MOND, Poisson's equation V 2 0n = A-wGp is replaced by 
the field equation 



/' 



IV0II 
a 



v<t> 



4%Gp, 



(1) 



where ||...|| is the standard Euclidean norm, and <f> is the 
gravitational potential for MOND. The gravitational accel- 
eration is g = — V</> just as the Newtonian acceleration 
gN = — V0N- For a system of finite mass, ||V0| — -> as 
||x|| — > oo, where x is the position vector relative to the sys- 
tem's centre of mass. The function n(y) is constrained by the 
theory only to the extent that it must run smoothly from 
niv) ~ 2/ at t/ <C 1 (the so-called "deep-MOND" regime) 
to /i(y) ~ 1 at j > 1 (the Newtonian regime), with the 
transition taking place at y « 1 (i.e., when ||V<^>|] is of order 
of the characteristic acceleration ao — 1.2 x lCP 10 ms -2 ). 
In spherical symmetry equation ^ reduces to Milgrom's 
(1983) original phenomenological relation 
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= gN, 
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in which g and gN are parallel, g ~ gN when g/ao 1 and 
g ~ y/aogN when g/a -C 1. 

Over more than two decades the theory has been quite 
successful, resisting several attempts at falsification (see, e.g. 
Sanders & McGaugh 2002; Bekenstein 2006; Milgrom 2008; 
Bekenstein 2009) . However, there are cases in which MOND 
appears to have difficulties in explaining the observed data: 
the X-ray and gravitational lensing properties of clusters of 
galaxies (e.g., The & White 1988; Sanders 2007; Natara- 
jan & Zhao 2008), and in particular the "Bullet" cluster 
(Clowe et al. 2006; Angus et al. 2007), the internal dynamics 
of dwarf spheroidal galaxies (Kleyna et al. 2001; Sanchez- 
Salcedo, Reyes-Iturbide & Hernandez 2006; Nipoti et al. 
2008; Angus 2008; Angus & Diaferio 2009), the phenomenon 
of galaxy merging (Nipoti, Londrillo & Ciotti 2007a), and 
the vertical kinematics of the Milky Way (Nipoti et al. 
2007b; Bienayme et al. 2009). Here we present a further test 
of MOND using the Milky Way Globular Clusters (GCs). 

Typical globular clusters are not ideal systems to test 
the MOND hypothesis because they are characterized by 
high stellar-mass surface density and internal acceleration 
larger than ao. However, there are several cases of less dense 
GCs, with internal acceleration comparable to or smaller 
than oo, which represent good candidates to test MOND 
(Baumgardt, Grebel & Kroupa 2005, hereafter BGK05). A 
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complication that arise in testing MOND with GCs is the so- 
called "external field effect" (Bekenstein & Milgrom 1984). 
In MOND the internal dynamics of a system is affected by 
the presence of even an uniform external gravitational field 
g cxt . As a consequence, the interpretation in the context of 
MOND of GC dynamics is typically difficult because one 
needs to account for the presence of the external gravita- 
tional field due to the Galaxy. An exception is represented 
by the GCs populating the outermost region of the Galactic 
halo, which are far enough to experience only a negligible 
acceleration from the Milky Way (</ e xt < ao), and therefore 
constitute an ideal laboratory to test MOND. 

The application of MOND to GCs has been studied 
by different authors by estimating the overall velocity dis- 
persions of GCs through the virial theorem (e.g. BGK05) 
or by deriving the their MOND velocity-dispersion profiles 
through the Jeans equations (e.g. Moffat & Toth 2008). As 
in Newtonian gravity, also in MOND the latter approach has 
the limitation that it is not based on self-consistent models, 
so that it is not guaranteed that there is a non-negative 
distribution function corresponding to the assumed density 
distribution. Haghi et al. (2009) solved this problem by ob- 
taining equilibrium MOND GC models as end-products of 
N-body simulations run with the numerical code n-mody 
(Londrillo & Nipoti 2008), based on the MOND potential 
solver described in Ciotti, Londrillo & Nipoti (2006). N- 
body models are useful when the external field is important 
and the GC is not spherically symmetric. When the external 
field is weak, the GC can be well represented by a spherical 
stellar system and it is possible to construct self-consistent 
models without resorting to simulations. 

In this paper we present self-consistent, spherically sym- 
metric models of stellar systems that are MOND steady- 
state solutions of the Fokker-Planck equation, which are the 
analogues of Newtonian King (1966) models. We use them 
to predict the MOND velocity-dispersion profiles of GCs be- 
longing to the outer Galactic halo, which can be used to dis- 
criminate between MOND and Newtonian gravity. In Sec- 
tion 2 we describe the theoretical basis of our models and 
discuss their structural and kinematic properties. Section 3 
is devoted to the application of the models to a sample of 
six GCs located in the external parts of the Galaxy. Finally, 
we summarize and discuss our results in Section 4. 



2 SELF-CONSISTENT SPHERICAL MODELS 
OF GLOBULAR CLUSTERS IN MOND 

To describe the phase-space distribution of a spherically 
symmetric star cluster in MOND we choose the steady- 
state solution of the Fokker-Planck equation proposed by 
King (1966), but with equation |T} replacing Poisson's equa- 
tion. The distribution function is 



f{r, v) = /o exp 



tjj-tpo 
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where the effective potential tp is the difference between the 
cluster potential at a given radius r and the potential out- 
side the cluster tidal extent tp = <f> — </> cx t, ipo is the cen- 
tral effective potential, /o is a scale factor, v c — ^/—2tp is 
the cluster escape velocity, and ctk is a normalization term 
which is proportional to the central velocity dispersion. The 



density p and the 3D velocity dispersion a v can be obtained 
by integrating the distribution function: 
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The above equations can be written in terms of dimension- 
less quantities by substituting 
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where po = p(0) is the central cluster density and 
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is the core radius (King 1966), so we get 
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The above expressions give the density and velocity- 
dispersion radial profiles as functions of the dimensionless 
potential W(f), which can be calculated by solving the mod- 
ified Poisson equation {1} with the boundary conditions at 
the centre 
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dr 
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In spherical symmetry equation |T} can be written in the 
form 
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or, using dimensionless quantities, 
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where £ = a^/aor c is a dimensionless parameter, which is 
smaller for systems closer to the deep-MOND regime. 

For a given choice of the (Wo,£) pair, W[f) has been 
obtained by numerically integrating equation (|13[l from the 
centre imposing the boundary conditions indicated in equa- 
tion (llip (see Appendix). The corresponding radial density 
and 3D velocity-dispersion profiles are then projected on the 
.jplane of the sky, giving the line-of-sight (LOS) velocity dis- 
persion ctlos as 



<?los(-R) : 
where 
E(R) = 2 
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is the surface density. The described procedure is straightfor- 
ward and produces self-consistent equilibrium models, trun- 
cated at a tidal radius n. 

A particular case of this family of models was studied by 
Brada & Milgrom (2000) , who constructed steady state King 
models of stellar systems in deep-MOND [fi(y) = y] regime 
as initial conditions for their N-body simulations of dwarf 
galaxies. Here we extend their models to general MOND 
cases in which the internal acceleration is not necessarily 
everywhere small as compared to ao. 



2.1 Limits of validity of the models 

The models presented in the above Section were constructed 
neglecting the effects of an external field on the modified 
Poisson equation (fTJ, so some considerations on the limits 
of validity of the models are necessary. The MOND field 
of a system with density distribution p(r) in the presence 
of the external field g cxt can be obtained by solving equa- 
tion {TJ with boundary conditions V</> — > g ex t for ||x|| — > oo 
(Bekenstein & Milgrom 1984). Computing such a field is in 
general a difficult task, because the presence of an exter- 
nal field breaks the symmetry of the system. However, the 
departure from spherical symmetry can be safely neglected 
in the external acceleration regime <? ex t <C a-o (Milgrom & 
Bekenstein 1987), so our spherically symmetric models are 
not valid when the external field is comparable to or larger 
than ciq. 

Another important issue is the concept of escape veloc- 
ity in MOND models. In the formulation presented in Sec- 
tion[2l the velocity distribution at a given radius is truncated 
at the local escape velocity v e = \f— 2$ (see also Sanchez- 
Salcedo & Hernandez 2007; Wu et al. 2008). This condition 
is no longer valid in the case of an isolated MOND model 
(<7cxt = 0), in which a particle with arbitrarily large speed 
remains bound to the cluster. In a more realistic case, when 
the system is surrounded by other objects (<? C xt > 0), the to- 
tal gravitational potential will admit a maximum value <j> ex t 
at a finite distance from the cluster centre (see Wu et al. 
2007). A cluster star located at a distance r from the cluster 
centre, with a speed v > v c , will reach the radius r t with 
non-null velocity, feel the external potential as dominant and 
never return. 

Summarizing, on the basis of the above considerations 
on the external field effect and on the escape velocity in 
MOND, our models are acceptable approximations of GCs 
when the external acceleration field lies in the range < 
g C xt -C a . 



2.2 Structural and kinematic properties of the 
models 

We describe here the structural and kinematic properties of 
the models presented in the previous Section. For a given 
interpolating function /i, a given choice of Wo and £ corre- 
sponds to a different model. Here (and below, when not spec- 
ified otherwise) we adopt the interpolating function (Famaey 
& Binney 2005) 



1 

i i i i | i i i i | i i i r 




log r (pc) 



Figure 1. Density (top panel), gravitational potential (middle 
panel) and internal acceleration modulus g (bottom panel) as 
functions of the distance from the cluster centre for a Newtonian 
King model with Wo = 6, log(M/MQ) = 5 and r c = 5 pc (solid 
lines) and of a MOND model with Wo = 9.5, log(M/M ) = 5 
and ^ = 0.53 (dashed lines). 
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Figure 2. Projected surface-density and LOS velocity-dispersion profiles of Wo = 6 MOND models with £ = 0.1 (dot-dashed lines), 
£ = 1 (dotted lines) and £ = 10 (dashed lines). For comparison, the profiles of Wo = 6 King Newtonian models are shown in both panels 
as solid lines. 
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Figure 3. Projected surface-density and LOS velocity-dispersion profiles of £ = 1 MOND models with Wo = 3 (dot-dashed lines), 
Wo = 6 (dotted lines) and Wo = 9 (dashed lines). The £ = 1 MOND isothermal sphere profiles (Milgrom 1984) are shown for comparison 
in both panels as solid lines. 



The knowledge of the cluster mass M allows us to determine 
r c , po and <tk using the relations 
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and 
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where 



r 2 p dr, 
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with f t = r t /r c . The above expressions can be used to trans- 
form the dimensionless quantities to physical ones. 
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In Fig. [T] we plot the density profile (top panel), the 
gravitational potential (central panel) and the acceleration 
modulus (bottom panel) of a Newtonian King model (solid 
curves) with Wo = 6, log(M/M ) = 5 and r c = 5 pc 
and, for comparison, the same quantities for a MOND model 
(dashed curves) having the same mass and reproducing the 
same surface-density profile (Wo = 9.5, f = 0.53). As can 
be seen, the MOND model has a potential well ~3 times 
deeper than the Newtonian model's one. Correspondingly, 
the MOND model has a substantially higher acceleration 
than the Newtonian model at all radii (see the bottom panel 
of Fig. [TJ. It must be noted that in the considered example, 
which is representative of a low surface-density cluster, the 
Newtonian acceleration is at all radii significantly lower than 
da, so the effects of MOND are large at any distance from 
the cluster centre. 

In Figs. [2] and [3] the surface-density and LOS velocity- 
dispersion profiles for a set of models with different choices of 
Wo and £ are shown. From Fig.[5]it is apparent that MOND 
models predict steeper density and velocity-dispersion pro- 
files with respect to Newtonian models having the same 
values of Wo- Indeed, the larger acceleration predicted by 
MOND at large distances from the cluster centre implies a 
steeper increase of the effective potential ip that reaches zero 
at a smaller radius (r t ) with respect to the Newtonian case. 
MOND models tend to Newtonian ones when large values 
of £ are considered (see Fig. [2J : the Newtonian case can be 
reproduced for ao — * (i.e. no characteristic acceleration). 
In this limit, £ — > oo, /J, — > 1, and equation (|13|l approaches 
Poisson's equation. Wo measures the depth of the gravita- 
tional potential well, so — like in Newtonian King models — 
as Wo —> oo, the escape velocity goes to infinity at any dis- 
tance from the cluster centre, the velocity distribution tends 
to the Maxwellian distribution, and the models approach the 
isothermal sphere (see Fig. [3}. 

We note that, with the exception of the isothermal 
sphere, in our MOND models ctlos goes to zero when r 
approaches the tidal radius (see right panels of Fig. [2] and 
Fig- EJ, as it must for the system to be bounded. In fact, as 
r — > rt the effective potential ip — > and the velocity distri- 
bution defined in equation © tends to a null width function. 
This is not the case for the isothermal sphere whose distribu- 
tion function has the same shape regardless of the potential, 
and it is not truncated, though having finite mass in MOND 
(Milgrom 1984). 



3 APPLICATION TO GLOBULAR CLUSTERS 
IN THE OUTER GALACTIC HALO 

The examples shown in the previous Section suggest that the 
shape of the velocity-dispersion profile of a GC of given den- 
sity profile can be quite different in MOND and Newtonian 
gravity. This makes the velocity-dispersion profiles of GCs 
a useful tool for distinguishing between MOND and New- 
tonian gravity. Here we apply our models to specific GCs, 
which may be used to test the MOND theory by combining 
measures of surface-brightness and LOS velocity-dispersion 
profiles. The models presented in this paper are valid in 
the external acceleration range < g cx t <C ao (see Sec- 
tion [271]) • In spite of this limitation, a number of GCs be- 
longing to the Galactic halo satisfy this constraint. Unfor- 



tunately, although the surface-brightness profiles of these 
clusters are well measured (McLaughlin & van der Marel 
2005, hereafter MvdM05), to date none of them have ac- 
curate velocity-dispersion profiles. A first attempt at esti- 
mating the velocity-dispersion profile of NGC2419 has been 
done recently by Baumgardt et al. (2009). 

Here we compare the velocity-dispersion profiles pre- 
dicted by the Newtonian and MOND models that repro- 
duce the surface-brightness profiles of six GCs located in 
the outer halo of the Milky Way (at distances > 50 kpc 
from the Galactic centre). These clusters (listed in Ta- 
ble 1) are subject to an external Galactic acceleration 
<?ext^3 x 10 -11 ms" 1 <C ao, so they are fully in the regime 
of validity of our models. 

3.1 MOND and Newtonian velocity-dispersion 
profiles for given mass and size 

Though the masses and physical sizes of GCs are known 
only with non-negligible uncertainties, it is useful to dis- 
cuss first the idealized case in which these quantities are 
given. In other words, for each object we fix a value of the 
mass-to- light ratio M/L and of the cluster distance, valid 
for both Newtonian and MOND models, which is used to 
convert surface brightness in surface mass density. We de- 
fer to Section T3.2I a discussion of the effects of varying M/L 
and distance. For Newtonian models we adopt the central 
dimensionless potentials Wo, core radii r c and masses M 
from MvdM05. For MOND models, we adopt the cluster 
mass estimated by MvdM05 and search for the values of Wo 
and £ that best reproduce the surface-density profile of the 
corresponding Newtonian models. 

In Table 1 the main structural properties derived for the 
six considered clusters are summarized. For Eridanus (which 
is not included in the catalog of MvdM05) we adopted the 
concentration, core radius, visual absolute magnitude given 
by Harris (1996) and estimated its mass by adopting a mass- 
V-band-luminosity ratio M/Lv = 1.892 (the same value 
adopted by MvdM05 for Pal 4, which has similar age and 
metallicity; Catelan 2000). In both Newtonian and MOND 
cases, the central LOS velocity dispersion ctlos.o has been 
calculated by using equation (|14|l . 

Figures [~4] [S] and show the surface- brightness profiles 
(left panels) and velocity-dispersion profiles (right panels) 
of the Newtonian and MOND models of the six consid- 
ered GCs. In the left panels we report also, where available, 
surface-brightness measures from Trager, King & Djorgov- 
ski (1995). As can be noted, MOND predictions significantly 
differ from Newtonian ones: in particular, it is apparent that 
the shape of MOND and Newtonian LOS velocity-dispersion 
profiles are remarkably different. As expected, MOND mod- 
els predict a larger velocity dispersion along the entire clus- 
ter extent with respect to Newtonian models with the same 
mass. The physical reason at the basis of this result is that, 
as showed in Fig. (TJ the same observed surface mass den- 
sity profile is fitted by MOND models with a deeper po- 
tential well and larger core radii. As a consequence, at a 
given radius, the corresponding escape velocity is larger and 
according to equation ((3]) the width of the velocity distri- 
bution turns out to be larger. The same qualitative results 
have been obtained by Moffat & Toth (2008) and Haghi et 
al. (2009), who adopted different methods to calculate the 
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Figure 4. Surface-brightness (left panel) and LOS velocity-dispersion (right panel) profiles according to MOND (dashed lines) and 
Newtonian (solid lines) models for the globular clusters NGC2419 and Eridanus. The surface-brightness measures for NGC2419 from 
Trager et al. (1995) are overplotted in the upper left panel with grey points. 



velocity-dispersion profiles of the sample of GCs indicated 
by BGK05. 

We note that the MOND ctlos profiles of our clusters 
are in general decreasing functions of radius and olos goes 
to zero when r approaches the truncation radius, consistent 
with the fact that the no stars can cross the system's bound- 
ary (see Sect. 12.2)) . An exception is the GC AMI, for which 
the predicted MOND velocity-dispersion profile in AMI is 
flat along the entire cluster extent. The surface-brightness 
profile of this cluster is indeed well fitted by the MOND 
isothermal sphere model (Milgrom 1984) corresponding to 
its mass. The maximum difference in the velocity dispersion 
predicted by MOND and Newtonian theories ranges from 
1.2 kms" 1 (Eridanus) to 2.4kms~ 1 (NGC2419), well above 
the accuracy currently achievable with high-resolution spec- 
troscopic analyses. 



3.2 Dependence on mass-to-light ratio and 
distance 

While the shape of the density and velocity-dispersion pro- 
files of Newtonian King models does not depend on the 
structural parameters (mass and core radius) , this is not the 
case for the corresponding MOND models. In these cases, 
the shape of the profiles varies for varying r c and/or M. 
In other words, while Newtonian King models of given Wo 
can be rescaled to represent systems of arbitrary values (in 
physical units) of r c and M, a MOND model of given Wo 
and £ represents only systems such that 

K =^L=^I{Wo,i) (21) 

(LOT c 

(see also Nipoti, Londrillo & Ciotti 2007c, for a detailed 
discussion of scaling of MOND models) . For each value of k 
there is a unique pair of parameters (Wo, £) that reproduces 
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Figure 5. Same as Fig. [4] but for the globular clusters AM 1 and Pal 3. 



Table 1. Structural parameters of the six considered outer halo globular clusters. 













Newtonian 






MOND 






Name 


log(M/M ) 


9cxt/ao 


M/Ly 


Wo 


r c 


CLOS.O 


W 




K 


°"LOS,0 


A max 








m q/ l v,q 




(PC) 


(km s- 1 ) 








(km s — 1 ) 


(km s" 1 ) 


NGC2419 


5.95 


0.12 


1.903 


6.5 


8.41 


5.31 


10.0 


1.25 


11.470 


6.59 


2.43 


Eridanus 


4.21 


0.11 


1.892 


5.3 


6.56 


0.97 


8.5 


0.10 


0.118 


2.13 


1.23 


AMI 


4.01 


0.09 


1.868 


6.6 


7.16 


0.61 


oo 


0.04 


0.028 


1.74 


1.69 


Pal 3 


4.65 


0.11 


1.869 


5.3 


10.68 


1.26 


8.5 


0.10 


0.118 


2.74 


1.58 


Pal 4 


4.58 


0.10 


1.892 


4.5 


12.16 


1.21 


7 


0.09 


0.079 


2.70 


1.51 


Pal 14 


4.09 


0.16 


1.885 


4.3 


19.60 


0.56 


6.5 


0.02 


0.004 


2.02 


1.46 



a given shape of the surface-brightness profile. The shape 
of the corresponding velocity-dispersion profiles is different 
for different values of re. This is illustrated in Fig. [7] where 
the velocity-dispersion profile of a Newtonian model with 
Wo = 6 is compared with a family of MOND models that 
share the same surface-brightness profile, but have different 



values of re. Here the LOS velocity dispersion and the radius 
are normalized to ctlos,o and r c , respectively, to highlight 
the different shape of the profiles. It is evident that while 
small values of re produce steep profiles, the MOND profiles 
approach the Newtonian one for increasing re. In fact, the 
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Figure 6. Same as Fig. [4] but for 

larger k the larger the internal acceleration of the cluster, 
which eventually exceeds ao over most of the cluster extent. 

In a practical application, once the surface-brightness 
profile of a given cluster is known, a unique pair of param- 
eters (Wo, that reproduces the shape of the observed 
surface-brightness profile can be determined only when an 
estimate of M and r c is provided. On the other hand, the 
uncertainties on the mass-luminosity ratio and distance pro- 
duce an uncertainty on M and r c and, consequently, on the 
predicted shape of the velocity-dispersion profile. 

To illustrate this issue, in Fig. [5] we show the velocity- 
dispersion profiles of the models that reproduce the surface- 
brightness profile of NGC2419 assuming a different mass 
(central panel) and a different core radius (top panel). In 
particular, we let the cluster mass vary by A log (M/Mq) = 
±0.15, with respect to the reference value log (M/Mq) = 
5.95 (therefore exploring the cases M/Lv ~ 2.7 and 
M/Lv = 1.35, beside the reference case M/Lv = 1-9), and 
we let the cluster core radius vary by A r c = ±0.1 r c (explor- 



i i i r 
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globular clusters Pal 4 and Pal 14. 

ing the cases of cluster distance d = 81 kpc and d — 99 kpc, 
beside the reference case d = 90 kpc). The corresponding 
overall cluster velocity dispersions (calculated by integrat- 
ing (Jlos over the entire cluster extent) are listed in Table 2. 
While a change in the core radius does not significantly af- 
fect either the shape or the central value of the velocity- 
dispersion profile, an even relatively small variation of M 
significantly alters the overall value of the velocity disper- 
sion. 

It must be noted that even comparing the highest M/L 
(and smallest distance) Newtonian model with the lowest 
M/L (and larger distance) MOND model, the resulting LOS 
velocity-dispersion profiles are remarkably different and thus 
observationally distinguishable. 



3.3 Dependence on the interpolating function fi 

An additional uncertainty is related to the adopted form of 
the interpolating function /j, (appearing in equation]!} , which 
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Figure 7. LOS velocity-dispersion profiles of a Newtonian King 
model with Wo = 6 (solid line) compared with a family of MOND 
models with different values of k, but with the same shape of the 
surface-brightness profile. 

is not constrained theoretically, except for the asymptotic 
behaviour (see Section [TJ. The choice of the functional form 
of n determines the behaviour of the MOND acceleration 
strength of the MOND effects in the intermediate accelera- 
tion regime (when g is of the order of ao), thus changing the 
overall shape of the velocity-dispersion profile. 

In the bottom panel of Fig. [8] we show the velocity- 
dispersion profiles of the models that reproduce the surface- 
brightness profile of NGC2419 assuming different forms of 
the interpolating function /j, (see Table 2), but keeping fixed 
M and r c . As expected, both the morphology of the velocity- 
dispersion profile and its average value depend on [i. How- 
ever, the effect of varying fx is typically small if we limit 
ourselves to standard proposal such as equation (|16[) (model 
MO) or Milgrom's (1983) n(y) = y/y/l + y 2 (model Ml), 
and we exclude unrealistic cases such as the step function 
adopted in model M2. 

3.4 Comparison with velocity-dispersion measures 

We have seen that, at least for the representative case of 
NGC2419, while the same overall velocity dispersion can 
be predicted by both Newtonian and MOND models with 
different choices of the mass-to-light ratio or of the interpo- 
lating function (compare, e.g., model Nm+ with models Ml, 
M2, and Mm- in Table 2), MOND profiles can be always eas- 
ily distinguished from Newtonian ones. The analysis of the 
shape of the velocity-dispersion profile represents therefore 
a robust method to discriminate between the two gravity 
theories. This method is indeed less sensitive to the errors 
on the cluster mass with respect to the simple comparison 
between the overall cluster velocity dispersion proposed by 
BGK05. Moreover, the approach suggested by these authors 
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Figure 8. LOS velocity-dispersion profiles for NGC2419. In the 
top panel Newtonian models (solid lines) are compared with 
MOND models (dashed lines) with different assumptions on the 
cluster distance d. In the central panel Newtonian models (solid 
lines) are compared with MOND models (dashed lines) with dif- 
ferent assumptions on the mass-to-light ratio. In the bottom panel 
the Newtonian model NO (solid line) is compared with MOND 
models MO (dashed line), Ml (dotted line) and M2 (dot-dashed 
line) differing in the choice of the interpolating function fi (see 
Table 2). The data by Baumgardt et al. (2009) are overplotted in 
the three panels. 
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Table 2. Predicted overall velocity dispersion for different 
MOND and Newtonian models of NGC2419. 



Model 


Kv) 


log (M/M Q ) 




< "LOS > 








(P c ) 


/ i ™ ~ — 1 \ 
( km s ) 


NO 


i 


5.95 


8.41 


4.40 


Nm- 


i 


5.80 


8.41 


3.70 


Nm+ 


i 


6.10 


8.41 


5.23 


Nr- 


i 


5.95 


7.57 


4.64 


Nr+ 


i 


5.95 


9.25 


4.20 


MO 


v/{i + v) 


5.95 


9.50 


6.18 


Ml 




5.95 


8.50 


5.56 


M2 


max(y, 1) 


5.95 


7.96 


5.22 


Mm- 


y/{i + y) 


5.80 


10.28 


5.51 


Mm+ 


v/0- + y) 


6.10 


9.23 


6.93 


Mr- 


V/Q- + V) 


5.95 


8.54 


6.28 


Mr+ 


y/(i + y) 


5.95 


11.08 


6.07 



favours low-mass GCs whose internal acceleration is lower 
than oo at any distance from the cluster centre. The low 
mass of these clusters, together with their large distances, 
imply a poor efficiency in measuring radial velocities for a 
meaningful sample of stars. 

Of course, measuring the velocity-dispersion profile is 
observationally more challenging than estimating the over- 
all velocity dispersion, but these kinds of measures are be- 
coming feasible even for relatively distant GCs. In Fig. [8] 
we overplot the three velocity-dispersion measures obtained 
by Baumgardt et al. (2009) from spectroscopic observations 
of 40 stars of NGC2419. Though there are large uncertain- 
ties due to the poor statistics, the overall trend defined by 
these observations appears hard to reconcile with MOND, 
at least under the considered assumption of spherical sym- 
metry and isotropic velocity distribution. This preliminary 
result confirms that NGC 2419 might be a crucial object to 
test MOND (as also suggested by Baumgardt et al. 2009) 
and strongly encourages future studies of this object com- 
bining higher resolution observations of a larger number of 
stars of NGC2419 as well as a systematic study of the possi- 
ble effects of orbital anisotropy, rotation and deviation from 
spherical symmetry. Among these effects, that of orbital 
anisotropy is likely the most important, because the shape of 
the LOS velocity-dispersion profile can depend significantly 
on the distribution of stellar orbits. In particular, a radi- 
ally anisotropic system is expected to have centrally steeper 
g los profile than an isotropic system with the same spatial 
distribution. Therefore, it is worth investigating whether ra- 
dially anisotropic MOND models can be reconciled with the 
velocity dispersion data for NGC2419. We address this ques- 
tion in the next Section. 



3.5 Radially anisotropic models of NGC2419 

Ideally, to explore the effect of orbital anisotropy on 
the olos profiles of GCs, one would need self-consistent 
anisotropic MOND models derived from the distribution 
function. Constructing such models is beyond the purpose 
of the present work: here we perform a preliminary analysis 
based on the numerical integration of the Jeans equations. 
We follow the standard procedure (Binney & Mamon 1982) , 
but with the MOND gravitational field replacing the New- 



tonian field. In practice, we take the spherically symmetric 
density distribution of the King model of NGC2419 (pa- 
rameters in Table 1), rescale it for the assumed value of 
M/Ly, and compute the MOND field generated by this 
density distribution using equation @. We then solve the 
Jeans equations assuming an anisotropy-parameter profile 
(3{r) = 1 — of (r)/2of(r), where of and of are the tangen- 
tial and radial components of the velocity dispersion ten- 
sor. Finally, we obtain the olosM by deprojecting oj?(r). 
For comparison, we also obtain olos profiles of anisotropic 
Newtonian models using the same procedure. 

As stressed in the Introduction, the Jeans-equations ap- 
proach does not guarantee that the obtained models are self- 
consistent. However, we can at least use some necessary con- 
ditions for consistency (e.g., Ciotti & Pellegrini 1992; An & 
Evans 2006; Ciotti & Morganti 2009) to exclude unphysical 
P{r). We note that these necessary conditions, though de- 
rived in the context of Newtonian gravity, apply to our self- 
gravitating MOND models, because each of this models can 
be formally interpreted as a non-self-gravitating distribution 
of tracer stars in the presence of a dominant mass distribu- 
tion having the same gravitational potential as the MOND 
potential of the cluster. An & Evans (2006) show that a 
necessary condition for consistencjQ is /3(0) ^ 7/2, where 
7 = — lim, — ,0 d hi p/dln r is the central logarithmic slope 
of the stellar density distribution. For NGC2419 7 ~ 0, so 
models with /3(0)>0 are inconsistent, implying that spheri- 
cal, radially anisotropic models with f3 independent of radius 
are unphysical. We then consider Osipkov-Merritt (here- 
after OM) models (Ospikov 1979; Merritt 1985), which are 
isotropic in the centre and radially anisotropic at large radii, 
having 

0(r) = p^*, (22) 

where r a is the anisotropy radius. A necessary condition for 
the consistency of OM models is that (r 2 + r't)p(r) is a non- 
increasing function of radius (Ciotti & Pellegrini 1992). In 
the considered model of NGC2419 this condition is satisfied 
for r a >9.43pc ~ 1.12r c , where r c is the core radius given in 
Table 1. In Fig.|9]we plot the olos profile for maximally ra- 
dially anisotropic (r a = 1.12r c ) MOND and Newtonian OM 
models of NGC2419 for the three reference values of M/Lv 
adopted in Fig. [8] as expected, radially anisotropic mod- 
els predict higher olos in the centre and lower olos in the 
outer regions than the corresponding isotropic models (mid- 
dle panel in Fig. [8} . While radially anisotropic Newtonian 
models reproduce well the observational data from Baum- 
gardt et al. (2009), MOND anisotropic models still tend to 
predict too high velocity dispersion. 

Given the large uncertainties in the observational 
data, it is not excluded that the velocity-dispersion pro- 
file of NGC2419 can be reproduced by a MOND radially 
anisotropic model with low stellar mass-to-light ratio. How- 
ever, it must be stressed that the models plotted in Fig.|9]are 
extreme cases: they satisfy only the necessary condition for 



1 This condition applies to a distribution of tracer stars in a 
potential well with finite depth (Evans, An & Walker 2009), as 
is the case for the MOND potentials of our models (see middle 
panel in Fig. [TJ. 
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Figure 9. Same as the middle panel of Fig. [8] but for radially 
anisotropic (Osipkov-Merritt) Newtonian (solid lines) and MOND 
(dashed lines) models of NGC2419. 



consistency, so they may have to be excluded as unphysical 
or unstable. 



4 SUMMARY AND CONCLUSIONS 

We have presented self-consistent dynamical models of stel- 
lar systems in MOND that can be used to study GCs in the 
outer Galactic regions. These models, which are the ana- 
logues of the Newtonian King models, are able to reproduce 
the observed surface-brightness profiles of GCs and can be 
used to predict their projected velocity-dispersion profiles in 
MOND. By comparing our models with the corresponding 
Newtonian ones, we found that it is impossible to reproduce 
simultaneously the same density and velocity-dispersion pro- 
files with both gravitational theories regardless of any choice 
of the free parameters of the models. This indicates GCs are 
one of the best laboratories to test the gravity theory in the 
low-acceleration regime, as already suggested by various au- 
thors (Scarpa, Marconi & Gilmozzi 2003; BGK05; Scarpa et 
al. 2007; Moffat & Toth 2008; Haghi et al. 2009; Lane et al. 
2009). 

Of course, the same effect produced by MOND on the 
velocity-dispersion profile can by reproduced by an ad hoc 
distribution of dark matter. The presence of dark-matter 
halos in GCs is predicted by some theories of GC forma- 
tion and evolution (see Mashchenko & Sills 2005 and refer- 
ences therein) and its observational evidence is still matter 
of debate (Moore 1996; Forbes et al. 2008). Thus, the de- 
tection of velocity-dispersion profiles in agreement with the 
predictions of MOND would not necessarily be a problem 
for the dark-matter paradigm. On the other hand, it must 
be stressed that the detection of velocity-dispersion profiles 



expected on the basis of Newtonian dynamics (without dark 
matter) could invalidate MOND. 

Testing MOND by using nearby GCs has been already 
attempted by Scarpa et al. (2003) who found that their 
velocity-dispersion profiles deviate from the prediction of 
Newtonian dynamics (without dark matter) at large dis- 
tance from their centres. However, the nearby clusters anal- 
ysed by these authors experience an external acceleration 
due to the Milky Way gravitational field that is larger than 
the critical acceleration a (BGK05; Moffat & Toth 2008). 
BGK05 indicated a sample of eight GCs where the predic- 
tions of MOND and Newtonian theories on the overall LOS 
velocity dispersion significantly differ. However, the method 
proposed by these authors is very sensitive to the adopted 
mass and distance of these clusters. Given the large dis- 
tances and low mass of these clusters, it is very difficult to 
observe the significant number of cluster member stars nec- 
essary to an accurate estimate of the velocity dispersion. In 
this regard, the difference between the MOND and Newto- 
nian predictions estimated by these authors never exceeds 
Aolos < 1.3kms _1 . Given the best accuracy achievable by 
the current observing facilities (~ 0.5 kms -1 ), it is hard to 
reach a firm conclusion on the validity of MOND using this 
approach. Nevertheless, the first results of the application of 
this method to Pal 14 suggests that the observed kinematic 
of this system might be a problem for MOND (Jordi et al. 
2009). 

A more robust test is to compare the observed shape of 
the velocity-dispersion profile with the prediction of New- 
tonian and MOND models (see also Moffat & Toth 2008; 
Haghi et al. 2009). We demonstrated that this method al- 
ways allows an unambiguous distinction between Newto- 
nian and MOND scenarios even when large uncertainties 
on the cluster mass and core radii are present. The best 
target GC, for which this approach is expected to be par- 
ticularly efficient, is NGC2419. Indeed, although not in- 
cluded in the list of BGK05, this cluster is the one with 
the largest absolute difference in the predicted velocity- 
dispersion profile by Newtonian and MOND models. More- 
over, it is not significantly affected by Galactic tidal ef- 
fects (Gnedin & Ostriker 1999) that can alter the shape 
of the velocity-dispersion profile in the outermost regions 
(Johnston, Sigurdsson & Hernquist 1999) and is massive 
enough to ensure a large number of stars at magnitudes 
easily reachable by the current generation of spectrographs. 
Recently Baumgardt et al. (2009) measured the velocity dis- 
persion at three different radii in this cluster. The two out- 
ermost data points are significantly lower than the velocity- 
dispersion profiles predicted by our isotropic MOND models 
of NGC2419, at least for stellar mass-to-light ratios in the 
range 1.35<Af/LvJs2.7. Only assuming quite strong radial 
orbital anisotropy and lower stellar mass-to-light ratio the 
velocity dispersion predicted by MOND can be reconciled 
with the observed data of NGC2419. Reproducing the ob- 
served kinematics of NGC2419 represents a challenge for 
MOND, although better data sets (larger number of stars 
and higher spectral resolution) and more sophisticated mod- 
elling are needed. 

A limitation of the present study is that our models 
are spherically symmetric, non rotating and with isotropic 
velocity distribution. The assumption of spherical symme- 
try and absence of significant rotation is in general justified 
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by the round appearance of GCs. In general, non-sphericity 
can be a problem in the determination of surface brightness 
and velocity dispersion profiles when large ellipticities are 
present (Perina et al. 2009). This effect should have only a 
minor impact in GCs which have generally symmetric den- 
sity contours, at least within few core radii. The amount of 
anisotropy in many GCs is estimated to be relatively small 
(Ashurov & Nuritdinov 2001), though a non- negligible de- 
gree of anisotropy is likely to be present in few GCs (Meylan 
& Heggie 1997 and references therein). Orbital anisotropy is 
also predicted by N-body simulations as a result of both 
primordial and evolutionary reasons (Giersz 2006). These 
effects are expected to be at least partially erased in Galac- 
tic GCs by the strong tidal interaction with the Milky Way 
which removes the initial velocity anisotropies and angular 
momentum making them more spherical (Goodwin 1997). 
A preliminary exploration of the effects of an extreme radial 
anisotropy shows a significant degeneracy between gravity 
law and orbital anisotropy. Nevertheless, the MOND and 
Newtonian models can be distinguished when high precision 
data are available. Another possible complication can be due 
to the presence of a significant fraction of unresolved binaries 
which can inflate the observed velocity dispersions (Cote et 
al. 2002). For instance, in NGC2419 the orbital velocity of 
a pair of equal-mass 0.8 Mq stars separated by a ^ 14 AU 
could be as high as 10 Km s _1 , remaining stable against 
collisional disruption (Hills 1984). A significant fraction of 
binaries can therefore increase the observed velocity disper- 
sion by few Km s _1 . This effect is particularly important 
in the central part of the cluster where binaries preferen- 
tially sink as a result of mass segregation. This could explain 
why the measured central velocity dispersion in NGC2419 
seems higher than the prediction of both Newtonian and 
MOND models. Note however that, in the case of NGC2419, 
the "binary-corrected" velocity dispersion would stray even 
more from the prediction of MOND models in the external 
region of the cluster. Given these uncertainties, a systematic 
exploration of the effects on the cluster kinematics of orbital 
anisotropy, rotation and deviation from spherical symmetry 
in general would be valuable also in MOND as well as in 
Newtonian dynamics (see Bertin & Varri 2008, and refer- 
ences therein). 
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APPENDIX A: NUMERICAL INTEGRATION 



We describe briefly the procedure used to numerically inte- 
grate equation (|13p . from f = to a given radius f , to obtain 
W(f). Equation (|13[) can be written as 



2 AW d 2 W 
f df df 2 



dW d 2 W 



df df 2 



dW 



df 



-9p(Al) 



where 



M'(y) ; 



dy 



(A2) 



The term d 2 W/df 2 diverges at the centre, but, for any 
choice of the MOND interpolating function fi, substitutions 
of power series show that 



lim 



d 2 W 



+o df 2 V 4£ 
whose integral converges and admits the exact solution 



(A3) 



3.i 



dW 

df 



wy^o = Wo - J^-f i. 

V 

We used the above equations to calculate the dimensionless 
potential and its derivatives at the centre. For the subse- 
quent integration steps we adopted a radial step Af such 
that 



( dW . _ d 2 W . _ . 
— Ar,— -Ar | < 0.01, 



(A4) 



in order to ensure a negligible numerical error in the inte- 
grations. 



This paper has been typeset from a TjtX/ P/IpX file prepared 
by the author. 



